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ABSTRACT 



By searching the IRAS and ISO databases we compiled a list of 60 debris disks 
^ ' which exhibit the highest fractional luminosity values {fd > 10~^) in the vicinity 

c3 I of the Sun {d < 120 pc). Eleven out of these 60 systems are new discoveries. 

Special care was taken to exclude bogus disks from the sample. We computed 
the fractional luminosity values using available IRAS, ISO, and Spitzer data, 
and analysed the galactic space velocities of the objects. The results revealed 
that stars with disks of high fractional luminosity often belong to young stellar 
kinematic groups, providing an opportunity to obtain improved age estimates 
for these systems. We found that practically all disks with > 5xl0~^ are 
younger than 100 Myr. The distribution of the disks in the fractional luminosity 
versus age diagram indicates that (1) the number of old systems with high fd 
is lower than was claimed before; (2) there exist many relatively young disks 
of moderate fractional luminosity; and (3) comparing the observations with a 
current theoretical model of debris disk evolution a general good agreement could 
be found. 



- 2 - 



Subject headings: circumstellar matter — infrared:stars — starsrkinematics 

1. Introduction 

One of the major discoveries of the IRAS mission was that main-sequence stars may 
exhibit excess emission at infrared (IR) wavelengths ("Vega-phenomenon", Aumann et al. 
1984). Systematic searches in the IRAS catalogues (Backman & Paresce 1993; Mannings & 
Barlow 1998; Silverstone 2000; Zuckerman & Song 2004a, and references therein) revealed 
that ~ 15% of main-sequence stars show infrared excess (Plets & Vynckicr 1999). It was 
suggested already after the first discovery that the excess can be attributed to thermal 
emission of dust confined into a circumstellar disk (Aumann et al. 1984). The existence 
of such debris disks was first confirmed by the coronographic observation of scattered light 
from the (3 Pic system (Smith & Terrile 1984). Subsequent imaging of specific systems at 
mid-infrared (e.g. HR 4796A, Koerner et al. 1998) and submillimeter wavelengths (e.g. e 
Eri, Holland et al. 1998) supported this picture. 

The possibility that debris disks might evolve over time was first mentioned by Backman 
& Gillett (1987) who proposed that disks around Vega, e Eridani, and Fomalhaut could 
be more evolved analogs of the /? Pic system. Using submillimeter measurements of the 
best known Vega-hke stars Holland et al. (1998) showed that the derived dust mass of the 
disks decay with stellar age as a power-law. Subsequent studies with the Infrared Space 
Observatory (ISO, Kessler et al. 1996) demonstrated that debris disks are more common 
around young stars (t < 400 Myr) than around old ones and there is a trend for older 
debris disks to be less massive than younger ones (Habing et al. 2001; Spangler et al. 2001; 
Silverstone 2000). The evolutionary picture was further refined by Decin et al. (2003), who 
reinvestigated the ISO results and revised the stellar age estimates. Recently Rieke et al. 
(2005) presented a survey of A-type stars, performed at 24/xm with the Multiband Imaging 
Photometer for Spitzer (MIPS, Rieke et al. 2004) onboard the Spitzer Space Telescope 
(Werner et al. 2004), and observed a general decay on a timescale of 150 Myr. 

A new generation of theoretical models has been developed to explain the temporal evo- 
lution of debris disks (Kenyon & Bromley 2002; Dominik & Decin 2003; Kenyon & Bromley 
2004) . These models take into account the fact that the destruction timescales of dust grains 
orbiting main-sequence stars are significantly shorter than the age of the central star, there- 
fore the observed dust grains in a debris disk must be continuously replenished (Backman & 
Paresce 1993). Collisional erosion of minor bodies in exosolar analogs of our Solar System or 
the sublimation of comets are the best explanations for the replenishment process (Harper 
et al. 1984; Backman & Paresce 1993). Most recent models link the temporal evolution of 
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debris disks to the formation and erosion of planetesimals (Kenyon & Bromley 2002, 2004). 

The amount of dust in a debris disk is usually characterised by its fractional luminosity, 
fd, defined as the ratio of integrated infrared excess of the disk to the bolometric luminosity 
of the star. Despite the general evolutionary trend described above, fractional luminosity 
values of individual systems were found to show a large spread of 10~^ ^ fd 10~^ at 
almost any age (Decin et al. 2003). Rieke et al. (2005) also found large variations of the IR 
excess around A-type stars within each age group. They emphasized the role of individual 
coUisional events between large planetesimals as one of the possible explanations of their 
result. 

Particularly interesting is the relatively high number of older systems (t>500 Myr) with 
high fractional luminosity values (/^ ^ 10^'^). These systems pose a serious challenge even 
to the new generation of theoretical models. A qualitative explanation for the existence 
of debris disks with large fa values around old stars was suggested by Dominik & Decin 
(2003). According to their theory, different planetesimal disks become active debris systems 
at different ages, because of the delayed onset of coUisional cascades. 

The presence of relatively old systems with high fractional luminosity, however, might 
partly be an observational artifact, since several factors could bias the distribution of debris 
disks on the fa vs. age diagram. Some stars have been nominated as Vega-candidates 
due to erroneous infrared photometric measurements, confusion by background sources, or 
the presence of extended nebulosity (where the IR emission is of interstellar rather than 
of circumstellar origin). In order to obtain a reliable picture of debris disk evolution these 
misidentifications ("bogus disks") have to be found and discarded. The actual positions 
of confirmed debris disks in the fd vs. age diagram might also be biased by measurement 
errors in the far-infrared photometry and even more significantly uncertainties in the age 
determination. 

In this paper we study a sample of debris disks which exhibit the highest fractional 
luminosity values in the solar neighbourhood. Setting a threshold value of fd = 10~^ and 
a distance limit of 120 pc we compiled a list of 60 disks, and performed an accurate deter- 
mination of their infrared excess using IRAS, ISO, and Spitzer data (Sect. 2). In Sect. 3 we 
present improved age estimates for a large number of stars, derived e.g., by determining their 
membership in young moving groups. In Sect. 4 we analyze the distribution of the disks in 
the fd vs. age diagram and find a significant fraction of our sample to be younger than 
was previously thought. Our conclusions are summarised in Sect. 5. In Appendix A. we list 
bogus disks identified in our work. Appendix B. lists new members of young moving groups 
discovered in the present study. 
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2. Sample selection 

We created the input list for this study by: (1) identifying debris disk candidates in 
the IRAS and ISO databases; (2) rejecting bogus debris disks and suspicious objects; (3) 
computing infrared fractional luminosity values and selecting disks with > 10"'*. 

2.1. Seeirching the IRAS catalogues for steirs with infrcired excess 

With the aim of compiling a list of main-sequence stars with IR excess, we made a 
systematic search in the IRAS Faint Source Survey Catalog (FSC, Moshir et al. 1989) 
and in the IRAS Serendipitous Survey Catalog (SSC, Kleinmann et al. 1986). In order to 
reduce source confusion, our survey was confined to \b\ > 10° Galactic latitudes. We selected 
all infrared sources with at least moderate flux quality at 25 or 60/xm, and their positions 
were correlated with entries from the Hipparcos Catalogue (ESA 1997) and the Tycho-2 
Spectral Type Catalogue (Wright et al. 2003). Positional coincidences within 30arcsec were 
extracted. In order to assure that the selected objects are not giant stars, the luminosity 
class was constrained to IV-V in the Hipparcos, and to V in the Tycho catalogue. We also 
included several objects whose luminosity class was not available in the Hipparcos catalogue 
but their absolute magnitudes indicated a main-sequence evolutionary phase. Since the 
infrared excess from early B-type stars might be due to free- free emission (Zuckerman 2001) 
the sample was limited to spectral types later than B9 in both catalogues. Our query is 
similar to that of Silverstone (2000) but - since we also considered the IRAS 25/xm band - 
it is also sensitive to stars with excess from warmer dust disk. 

Following the principles of the method by Plets & Vynckier (1999), for each selected 
star we predicted the far-infrared flux density of the stellar photosphere using the Kg-band 
magnitude (or V-band, when good quality Kg photometry was not available) and the B-V 
color index. Kg-band photometry was drawn from the Two Micron All Sky Survey (2MASS) 
catalog (Cutri et al. 2003), V magnitudes and B-V color indices were taken from the Hip- 
parcos and Tycho Catalogues. As a first step a photospheric 25//m flux density was derived 
from the Kg magnitude and the B-V color of the star using the collection of stellar model 
predictions by M. Cohen and P. Hammersley (available on the ISO Data Centre home page). 
Then, color relationships predicting the photospheric flux ratios between 25yLim and a selec- 
tion of IRAS, ISO, and Spitzer photometric bands were also derived from the same stellar 
models. The average accuracy of the predicted far-infrared fluxes is estimated to be around 
4% when computed from the X^-magnitudes, and 8% when computed from V-magnitudes. 

In order to compute IR excess values the predicted photospheric flux densities were 
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subtracted from the measured flux densities in each IRAS band. In principle the IRAS 
fluxes have to be color corrected since the shape of the spectral energy distribution of the 
system usually differs from the Fj,~z/~^ reference spectrum (this spectral shape was assumed 
while the flux densities quoted in the IRAS catalogues were derived from the detector in-band 
powers). Since the true spectrum of the system is not known a priori, we decided to multiply 
the predicted photospheric fluxes - rather than dividing the IRAS flux densities - with color 
correction factors appropriate for a stellar photosphere (IRAS Explanatory Supplement, 
Beichman et al. 1988). The significance level of the infrared excess was calculated in each 
photometric band with the following formula: 

o ^meas Fpred \ 

^excess i \^) 

where SF^ngas is the quoted uncertainty in the FSC or SSC, and dF^^ed is the uncertainty of 
the prediction described above. When Sexcess was greater than 3 either in the 25 or 60/xm 
bands, the object was selected as a excess candidate star. Applying the above criteria we 
identified in total 355 excess candidate stars in the IRAS databases. 



2.2. ISO-based selection of stars with infrared excesses 

In a second step the IRAS-based hst was supplemented with excess stars selected from 
the ISO databases. The Vega-phenomenon was a key programme for ISO (see Sect. 1) and 
a number of stars have been observed with ISOPHOT, the onboard photometer (Lemke et 
al. 1996). We collected all ISOPHOT observations of normal stars from different observers 
performed in mini-map, sparse-map or staring mode (a detailed description of these observing 
modes is given in The ISO Handbook Vol. IV, Laureijs et al. 2003), and performed a 
homogeneous re-evalution of the whole sample (Abraham et al. 2003; Abraham et al., 
2006, in prep.). For details of the data analysis and criteria for candidate excess stars, see 
Abraham et al. (2006, in prep.). We note that most selected candidates have already been 
pubhshed by the original observers (Decin et al. 2000; Habing et al. 2001; Spangler et al. 
2001; Silverstone 2000), but due to their different processing schemes the published flux 
densities cannot be directly merged for a homogeneous catalogue. The merged IRAS- and 
ISO-based lists includes altogether 364 IR excess stars. 
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2.3. Rejection of suspicious objects 

Since our goal is to compile a list of debris disks, we excluded all known young stellar 
objects (e.g. T Tauri or Herbig Ae/Be stars) which harbour protoplanetary disks. The 
sample could also be contaminated by source confusion: due to the low spatial resolution of 
IRAS at far-IR wavelengths, many of the positional coincidences between a star and a far-IR 
source could be bogus and the far-IR emission is related to a foreground or background 
object. 

For part of the sample (110 stars) higher spatial resolution infrared maps are available, 
obtained either by the ISOPHOT or the MIPS instrument. We downloaded ISOPHOT data 
from the ISO Data Archive (IDA) and processed with the Phot Interactive Analysis (PIA) 
version 10.0 (Gabriel et al. 1997). MIPS Basic Calibrated Data (BCD) files were downloaded 
from the Spitzer Science Center (SSC) data archive. These latter products are composed of 
2-dimensional FITS image files which included all general calibrations and corrections for 
MIPS detectors (Gordon et al. 2005). In each case these data were coadded and corrected for 
array distortions with the SSC MOPEX (MOsaicking and Point source Extraction, Makovoz 
& Marleau 2005) software. Bad data flagged in the BCD mask flies, as well as permanently 
damaged pixels flagged in the static pixel mask file were ignored during the data combination. 
Output mosaics had pixels with size of 2.5" at 24/im and 4" at 70//m. SSC MOPEX/APEX 
software package was used to detect sources and determine their positions on the final maps. 

The positions of the infrared sources were determined on the ISO and Spitzer maps and 
objects whose coordinates differed from the optical position (and in some cases coincided with 
a nearby background object), and/or associated with extended nebulosity, were discarded 
from the fist. For positional discrepancy the threshold value was set to half of the width of a 
point source's footprint. In the case of ISOPHOT the large pixel size dominated the footprint 
and in the 60-100/im range the threshold was 23". In the case of MIPS arrays footprint was 
defined by the telescope Point Spread Function and we adopted threshold values of 3" and 9" 
at 24 and 70fim, respectively. The absolute pointing uncertainty was less than these values 
for both satellites. In total 24 disk candidates were dropped from the list. 

When neither ISO nor Spitzer maps were available we made an attempt to filter out 
bogus disks by assuming that an object is possibly affected by source confusion if: 

• a known galaxy or evolved star (OH/IR source, Mira variable) is located within 1 
arcmin of the IRAS position; 

• a source included in the IRAS Small Scale Structure Catalog (Helou & Walker 1988) 
or in the 2MASS Extended Source Catalog (Jarrett et al. 2000) is located within 1 
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arcmin of the IRAS position; 

• a 2MASS source with an excess in the Kg band (identified in the H — Kg vs. J — H 
diagram in comparison with the locus of the main-sequence and taken into account the 
reddening path) is located within 1 arcmin of the IRAS position; 

• the 60-to-lOO/im flux ratio of the candidate source resembles the color of infrared 
cirrus < 0.25, which correspond blackbody temperatures lower than 33 K ). At 
least moderate flux quality flags were required in both IRAS bands. 

These cases were also discarded from our hst of debris disk candidates (48 objects). 

In an earlier study, Kalas et al. (2002) used high-angular resolution coronographic ob- 
servations at optical wavelengths and found cases where the far-infrared excess observed 
by IRAS was of interstellar - rather than circumstellar - origin ("Pleiades-phenomenon"), 
leading to false entries in the Vega-candidate lists. They also suggested that a signiflcant 
fraction of Vega-candidates beyond the Local Bubble might be bogus, since the star illumi- 
nates nearby interstellar matter rather than a circumstellar disk. After checking the positions 
of our sources projected on recent maps of the Local Bubble (Lallcment et al. 2003), we dis- 
carded all objects situated in the wall of the bubble or beyond. The wall was defined as 
the isocontours corresponding to the 50 mA Nal D2-line equivalent widths in the maps. In 
practice nearly all of our sources beyond 120 pc were removed, while within this radius only 
a few were dropped. Thus we defined a maximum distance limit of 120 pc for our stars, 
constructing a nearly complete volume-limited sample. 

2.4. The list of disks with high fractional luminosity 

In order to compute fractional luminosity values for each candidate star, we constructed 
spectral energy distributions by combining infrared fluxes from the FSC, SSC, ISOPHOT (re- 
evaluated by us, see Abraham et al., 2006, in prep.) and additional MIPS and submillimeter 
fluxes from the hterature. The excess above the predicted photosphere was fitted by a 
single temperature modified blackbody, where the emissivity was assumed to vary as 1 — 
exp [— (Ao/A)^] where Aq was set to 100/xm (see e.g. WiUiams et al., 2004). We fixed P 
equal to 1, which is a typical value in the case of debris systems (Dent et al. 2000). If 
the excess was detected at one wavelength only, we adopted a modified blackbody whose 
peak (in F,y) coincided with that single wavelength. From the fitted spectral shape color 
correction factors were computed and applied to the data. Then again a modified blackbody 
was fitted resulting in new color correction factors, and this procedure was repeated until 
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the color correction factors converged. Finally, the fractional dust luminosity was calculated 
as fd — Lm/Lhoi- In order to estimate the uncertainties on our fractional luminosity values 
we performed a Monte Carlo simulation. We added Gaussian noise to the photometric data 
points using their quoted la photometric errors and then recomputed the fractional dust 
luminosities. Formal uncertainties of the predicted theoretical photospheric fluxes were also 
taken into account. Final uncertainties were derived as the standard deviation of these 
values after 1000 repetitions. We note that these values include only random uncertainties; 
systematic errors due to e.g. limited wavelength coverage are not taken into account. 

Artymowicz (1996) argued that debris disks are confined to fd < 10^^ and sources with 
higher fractional luminosity probably contain a significant amount of gas (e.g. T Tau and 
Herbig Ae/Be stars, "transition" objects). Therefore we excluded objects with > 10~^ 
from our sample. Then the remaining sample was sorted by decreasing fa values and stars 
with fd > 10~^ (60 stars) were taken for the further analysis presented in this paper. 

Basic stellar parameters for these 60 objects, as well as derived fractional luminosities 
and their uncertainties, are presented in Table 1. Infrared data used in our analysis, includ- 
ing both the original flux values as listed in the catalogues or provided by our reduction 
algorithm, as well as corrected fluxes, where color correction was applied, are given in Ta- 
ble 2. The table also contains photospheric flux predictions for the specific wavelengths. 
Inspecting the flux density values obtained by different instruments at the same wavelength 
(e.g. IRAS and ISOPHOT at 60/xm) one finds discrepancies which may arise e.g. from the 
different beams and different calibration strategies of the instruments. Comparing the IRAS 
and ISOPHOT fiux values in Table 2 a general good agreement within la was found, with 
no deviations above the 3a limit (the uncertainty a was computed as the quadratic sum 
of quoted uncertainties from the two instruments). We also compared our ISOPHOT fiux 
densities with the results of earlier evaluations of the same observations in the literature. In 
most cases the results were consistent (except HD 10647 and HD 53143 where Decin et al. 
(2000) derived significantly higher values at 60/um; the probable explanation is that we used 
a more advanced version, VI 0.0, of the PI A software). 

As was discussed in Sect. 2.1, during the analysis we rejected several systems as bogus 
disks or suspicious objects. Those rejected stars which were previously proposed to harbour 
debris disks in the literature and would have been included in our final list (on the basis 
of their quoted fractional luminosity in the original paper) are presented in Appendix A. 
together with a brief description of the reason of rejection. 
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3. Age determination 

3.1. Membership in young moving groups 

Age determination for main-sequence field stars is cfiallenging, and sometimes results 
in very uncertain values. Ages of open cluster members, however, can be estimated more 
accurately e.g. by fitting their main-sequence locus in the color-magnitude diagram with 
theoretical isochrones or by determining the location of the "lithium depletion edge" in the 
cluster and comparing it with the predictions of theoretical evolutionary models. A number 
of young clusters (a Per, Pleiades, Hyades etc.) have been dated so far (Mcynct ct al. 1993; 
StaufTcr ct al. 1998). Similarly, the ages of young stellar kinematic groups, discovered mainly 
in recent years, are relatively well determined (e.g. Zuckerman & Song 2004b). It was a 
very important result that several stars with the strongest infrared excess turned out to 
be members of such moving groups, and in some cases the ages of these stars had to be 
revised significantly (e.g. the case of /3 Pic, Barrado y Navascues et al. 1999). In order 
to obtain more reliable ages for our sample, we performed a systematic investigation of the 
possible relationship between our excess stars and nearby kinematic moving groups, stellar 
associations or open clusters. 

A common method to decide whether an object belongs to a moving group is to compare 
its Galactic space velocity components with the mean velocity components of the group. 
In order to compute the space velocity for the stars in Table 1, we collected parallaxes 
and proper motions from the Hipparcos and Tycho-2 catalogues. When accurate parallax 
information was not available, a photometric distance was adopted. 

Radial velocities were taken from the literature (see Table 3 for references), or from our 
own observations. The new observations were carried out with the 2.3m ANU-telescope at 
the Siding Spring Observatory, Australia, on 11 nights between 21 March and 22 August 
2005. The spectra were taken with the Double Beam Spectrograph using a 1200mm~^ 
grating in the red arm. The recorded spectra covered 1000 A between 5800-6800 A, with a 
dispersion of 0.55 Apx~^. This leads to a nominal resolution of about 1 A. The exposure time 
ranged between 30 s and 200 s depending on the brightness of the target and the weather 
conditions. We obtained on average 4-7 spectra for each star. Since all the target stars are 
bright objects (V<10 mag), we could easily reach S/N~150-200 for every spectrum. All 
spectra were reduced with standard tasks in IRAF^. Reduction consisted of bias and flat 



-'^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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field corrections, aperture extraction, wavelength calibration and continuum normalization. 
We did not attempt flux calibration because the conditions were often non-photometric 
and the main aim was to measure radial velocities. Radial velocities were determined by 
cross-correlation, using the IRAF task fxcor, choosing HD 187691 as a stable lAU velocity 
standard. The cross-correlated region was 100 A centered on the Ha line, which is far the 
strongest spectral feature in our range. The finally adopted velocities were calculated as 
simple mean values of the individual measurements. Our experiences have shown that the 
typical measurement errors were about 4-7 km s~^ per point, so that the mean values have 
±1-3 km standard deviations. These were adopted as the uncertainties shown in Table 3. 

In the calculation of the Galactic space velocity we used a right-handed coordinate 
system (U is positive towards the Galactic centre, V is positive in the direction of galactic 
rotation and W is positive towards the North galactic pole) and followed the general recipe 
described in "The Hipparcos and Tycho Gatalogues" (ESA 1997). The computed Galactic 
space velocity components and their uncertainties are given in Table 3. 

Table 4 summarizes the basic properties of the relevant moving groups and associations 
within 120 pc from the Sun. The probability that star i is a member of moving group j, can 
be computed by: 



Pij = exp 



H — h 



(2) 



where Ui,Vi, Wi and o";/., cry., are the Galactic space velocity components of the star and 
their uncertainties, respectively, while Uj,Vj,Wj and auj,avj,crwj are the mean Galactic 
space velocity components of the specific kinematic group as well as the corresponding errors. 
In this formula we assumed that the velocity distribution within a group is Gaussian. When 
, , <ywj parameters were not available in the literature, we computed them from the 
velocity dispersion of known members around the mean. Our newly calculated mean values 
were always consistent with those from the literature within the uncertainties. In those 
few cases where no sufficient membership information could be found in the literature, we 
adopted au. — uy. — aw, — 2km s~^ (a characteristic value in the previous cases). 

Probability values for each star with respect to each group were computed. Then we 
checked the resulting values for objects already assigned to a group in the hterature. The 
numbers spread in the range of 0.2 < P < 1.0, therefore we set P = 0.2 as a lower limit 
for the new moving group member candidates as well. Stars assigned to any group above 
this threshold were further checked by comparing their 3-dimensional space location with 
the volume occupied by the group (most groups are rather confined in space). There were a 
few stars which could be assigned to both the Tucana-Horologium and GAYA2 associations; 
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these cases are analyzed in Appendix B. Table 3 presents the final assignments between stars 
and kinematic groups. 

Prom our sample of 60 objects, 26 sources could be linked to stellar associations; 13 
of them are new members identified in the present study. For 10 stars out of these 13, 
age estimates are available in the literature. For these ten objects we directly compare our 
age estimates with previous values in Table 5. In most cases ages derived in our study are 
younger than the earlier values. The discrepancy is particularly obvious in the case of ages 
derived from isochronc fitting (for example Nordstrom et al. 2004). However, this problem 
is not related only to the present study. For example, HD 105 has a moving group age 
of 30l^°Myr (Mamajek et al. 2004a) while Nordstrdm ct al. (2004) quoted 8600+^°°° Myr 
from isochrone fitting. Similar objects in our sample are HD 25457 (50 — 100 Myr from 
moving group, Zuckerman & Song 2004b and 4000l2ioo^y^ from isochrones, Nordstrom 
et al. 2004); HD 164249 (12l^Myr, Song et al. 2003 and 22001}^^^ Myr, Nordstrom et al. 
2004); HD 181327 ( 12l^Myr, Song ct al. 2003 and 1300t Myr, Nordstrom et al. 2004). 
The age uncertainty related to isochrone fitting might arise from the lack of information on 
whether the star is in the pre-main-sequence phase of its evolution or is an evolved object 
above the main-sequence. In ambiguous cases we always adopted age estimates derived from 
stellar kinematic group membership. 

3.2. Statistical age estimates for the disk sample 

For stars not assigned to any moving groups other age estimation methods are needed. 
Before focusing on individual systems, in this subsection we analyze what can be learned 
about the age distribution of our sample of debris disk systems. 

3.2.1. Distribution of the excess stars in the velocity space 

The distribution of the derived Galactic space velocities (Table 3) are displayed in 
Fig. la-b. Overplotted is the box occupied by young disk population stars defined by (Leggett 
1992) on the basis of a systematic study of Eggen (1989). The plots show that most stars 
from our sample belong to this population. This fact suggest that the majority of our sample 
of stars from Table 1 are relatively young. 
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3.2.2. Location of A-type stars on the CMD 

For a sample of bright A-type stars Jura et al. (1998) demonstrated that three objects 
with strong infrared excess (HR4796A, HD 9672, j3 Pic) are located close to the lower bound- 
ary of the distribution in the Color-Magnitude Diagram (CMD). Comparing their locations 
with the loci of the youngest, nearby open clusters (a Per, IC2391, Pleiades), Lowrance et al. 
(2000) argued that stars with strong infrared excess are typically younger than these clusters. 
Adopting this idea we selected stars with B— V ranging between -0.1 and 0.33 (corresponding 
mainly to A-type stars) from Table 1, and plotted them in the CMD of Fig. 2. In addition 
we overplotted a volume-limited sample of A-type stars {d < 100 pc) extracted from the 
Hipparcos Catalogue (it was requested that the parallax error was less than 10% and the 
B— V uncertainty was lower than 0.01 mag). For comparison, members of the open clusters 
a Per (80 Myr) and Hyades (600 Myr) are marked. The older cluster, Hyades, covers the 
upper part of the distribution, while a Per stars are situated at lower absolute magnitudes 
for the same color. Pleiades (100 Myr, not plotted in the figure for clarity) occupies the same 
region as a Per (see Fig. 3 in Lowrance et al. 2000). 

Figure 2 shows that the majority of the objects from our sample of high stars appear 
to be close to the lower boundary of the area occupied by A-stars, and are located below the 
region of a Per and Pleiades, with some overlap. This suggests that the early spectral type 
stars from our sample of high disks are young, probably close to the zero age main-sequence 
(ZAMS), and they are presumably not older than 100 Myr, the age of the Pleiades. 

3.3. Ages of individual objects 

For those objects which could be assigned to one of the moving groups or associations the 
age of that group as well as its uncertainty was adopted (26 objects). For a number of stars 
not associated with groups or with associations, age estimates could be found in the literature 
(19 stars, for references see Table 1). Sometimes literature data for a specific star scatter 
significantly; in these cases we adopted an age range which covers all quoted values and their 
uncertainties. When the literature search did not yield any dating, we made age estimates 
by plotting the stars on the HR diagram and comparing their positions to isochrones. For 
stars with spectral types in the range B9-G5 the isochrone age was estimated following the 
general outline described by Lachaume et al. (1999), using the Padova theoretical isochrones 
(Girardi et al. 2000). This method was applied to 13 objects. 

This isochrone method gives only upper limits for some A-type stars. As a best estimate 
for these stars (five cases), we adopted an upper limit of 100 Myr, consistent with our results 
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shown in Fig. 2 and discussed in Section 3.2.2. For stars of later spectral type there 
are some widely used age indicators, like the strength of the CallH&K lines or the X-ray 
luminosity of the star. In the case of HD 121812 our age estimate is based on the former 
method, taking the measured value from Strassmeier et al. (2000) and using the calibration 
of Lachaume et al. (1999). HD 130693 has a ROSAT counterpart and its X-ray luminosity of 
logLx = 29.7 erg s^^ was compared with the X-ray luminosity distribution function of late- 
type members in different associations (see figure 2 in Stelzer & Neuhauser 2000), yielding 
an age range of 10 — 100 Myr. In Table 1 we summarize the age estimates for each object. 



4. Discussion 

4.1. Connection between debris disks and young moving groups 

From our sample of 60 main-sequence stars exhibiting strong infrared excess, 26 can be 
assigned to young stellar kinematic groups. In order to test whether the frequency of stars 
belonging to young stellar kinematic groups is similar in a general sample we determined the 
corresponding ratio within a volume-limited sample of normal stars. First we created this 
sample by selecting stars from the Hipparcos catalogue using the following criteria: (1) they 
arc closer as 120 pc (the same volume-limit than in our sample); (2) their Survey flag in the 
catalogue^ was set to "S"; (3) they have radial velocity measurement with uncertainty less 
than 5 km s~^. The second condition guarantees that stars observed in various individual 
projects do not introduce a bias in the analysis of the velocity distribution (Skuljan et 
al. 1999). In addition, Binney et al. (1997) noted that radial velocities are preferentially 
observed for high-proper-motion stars which may cause a kinematical bias in our sample. 
Following the proposal of Skuljan et al. (1999) we constrained the sample for stars exhibiting 
low transverse velocity, and excluded all stars with Vt > 80 km s~^ in order to avoid this 
bias. 

The query resulted in 7519 objects, for which we computed the UVW Galactic space 
motion components. For each star in the two samples (the 60 debris systems in Table 1, and 
the newly deflned volume-limited stellar sample) we determined the Euclidean distance in 
the 3D velocity space from the closest moving group, dmin- In this analysis we considered 
only groups younger than 150 Myr. In Fig. 3 we plotted the histograms of 6min for the 



^Field H68 in the Hipparcos catalogue. 'S' indicates that the entry is contained within the 'survey', which 
was the basic Ust of bright stars added to and merged with the total list of proposed stars, to provide a 
stellar sample almost complete to well-defined limits. The limiting magnitude was a function of the stars's 
spectral types and galactic latitude. 
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two samples. A two-sided Kolmogorov-Smirnov test shows that the two distributions are 
different with a probabihty higher than 99.99%. This result indicates that debris systems 
of high infrared fractional luminosity are much more intimately linked to the nearby young 
stellar kinematic groups than the majority of normal stars. 

4.2. The relationship between fractional luminosity and age 

Zuckerman & Song (2004b) hypothesised that stars with > 10~^ are younger than 
100 Myr, and therefore a high fd value can be used ge indicator. This proposal is in 

contradiction with the conclusion of Dccin ct al. (2003), who claimed the existence of high 
fd disks around older stars. In order to test which proposal is supported by our data, in 
Fig. 4 wc plotted the distribution of ages as a function of the fractional luminosity fd from 
Table 1. We plotted in red the debris disks whose presence was explicitely confirmed by an 
instrument independent of IRAS (see Col. 8 in Table 1). The confirmation could be based 
e.g. on high spatial resolution infrared images (ISO/ISOPHOT, Spitzer/MIPS), or on mid- 
infrared spectra (Spitzer/IRS), or on coronographic images (HST/ACS, HST/NICMOS). We 
found 43 confirmed debris disks in total. 

Most data points with fd>5x 10~^ fall below the age threshold of 100 Myr (marked 
by a dashed line), while objects with lower fd show a larger spread in age. This trend 
can be recognized in the whole sample, but is especially clear in the confirmed subsample 
(red symbols). There is only one noteworthy case: HD 121812 with an age of 230lgo°Myr 
exhibits fractional luminosity exceeding the 5 x 10~"^ threshold value. However, the presence 
of a debris disk around this star has not been confirmed independently of IRAS. On the 
basis of this result we conclude that - according to the suggestion of Zuckerman & Song 
(2004b) - the majority of debris disks with > 5 x 10"*^ are younger than 100 Myr, and 
high fractional luminosities can be used as an indicator of youth. Nevertheless the opposite 
is not true, i.e. a low fd value is not correlated with age, and in particular is not an indicator 
of antiquity. 

There is a growing list of debris disk systems which have been discovered by the sensitive 
detectors of the Spitzer Space Telescope (Beichman et al. 2005a; Bryden et al. 2005; Chen et 
al. 2005a,b; Kim et al. 2005; Low et al. 2005; Meyer et al. 2004; Stauffer et al. 2005; Uzpen 
et al. 2005), and one may wonder whether these new observations support our previous 
conclusion. As a preliminary check we collected from the cited papers all debris disks with 
fd > 5x 10~^. We used fractional luminosity values and age estimates as quoted in the papers. 
We found that all of these new disks discovered so far belong to the ~16Myr old Lower 
Centaurus Crux subgroup of the Scorpius-Centaurus association (HD 106906, HD 113556, 
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HD 113766, HD 114082, HD 115600, HD 117214; Chen et al. 2005a), or to the TWA (TWA 7, 
TWA13A, TWA13B; Low et al. 2005), or to the star forming region RCW49 (18 possible 
warm debris disks, Uzpcn ct al. 2005), which unambigously shows that these objects are 
young, in agreement with the conclusion of the present paper. 



4.3. Debris disk evolution and the cases of old systems 

There are a number of models in the literature (see Sect.l) to describe the temporal 
evolution of debris disks. In the following we compare our results with predictions. 

Dominik & Decin (2003) proposed a simple coUisional model which assumes that all 
dust grains in the debris disk are produced in collisions between planetesimals within a ring 
whose radius is constant during the whole evolution. In coUisional equilibrium - when the 
dust production and destruction rates are in balance - the grain loss mechanism governs 
the amount of dust visible in the system. If dust destruction is dominated by collisions 
the fractional luminosity fd decreases proportionally to t'^. If the dust removal process is 
dominated by the Poynting-Robertson drag, fd oc It is predicted that in disks with 

fd > 10^'^ the evolution is dominated by collisions (Dominik & Decin 2003; Wyatt 2005). 

Three families of disk evolution models, computed from eqs. 7, 35-40 of Dominik & 
Decin (2003), are plotted as shadowed bands in Figure 5. The data points and symbols in 
this figure are identical to those in Fig. 4. The main differences between the model families 
are related to disk mass: (a) Ma = 10 M®; (b) Ma = 50^1®; and (c) Ma =250 M®. The 
width of each band corresponds to a range in stellar mass from 0.5Mq to 3Mq. Additional 
parameters are: the characteristic radius of the ring of planetesimals Tc = 43 AU; the radius of 
planetesimals Oc = 10 km; the density of the planetesimal material pc — 1.5gcm~^ (proposed 
for icy bodies with small rocky component, Greenberg 1998; Kenyon 2002); the size of the 
smallest visible grains a^ig = lOfim. (taken from Jura et al. 2004); the absorption efficiency 
of the dust particles Qabs = 1; and eg = 226 (defined in eq. 22 in Dominik & Decin 2003). 

Figure 5 shows that the location of most stars on the evolutionary diagram can be 
explained by the models. The number of older stars exhibiting high fd values incompatible 
with the models is relatively low. Most of these disks are located in the t > 10^ yr and 
/rf ^ 5 X 10^'^ area. A possible explanation for the origin of these stars was proposed 
by Dominik & Decin (2003) who assumed that different planetesimal disks become active 
debris system at different ages because of the delayed onset of the coUisional cascade. Large 
coUisional events may also increase the brightness of a debris disk temporarily (Rieke et al. 
2005). Extraordinary events during the evolution like e.g. a proposed 'supercomet' in the 
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HD 69830 system (Beichman et al. 2005b) cannot be excluded, too. Nevertheless, the low- 
number of systems incompatible with the models, especially at > 5 x 10~^, indicates that 
the above scenarios do not represent the main evolutionary trend. 

It is important to note that a large spread in fractional luminosities (10~^ < fa < 
5 X 10~^) can be observed in the figure among young debris systems {t < 100 Myr) . This result 
resembles the findings of Ricke et al. (2005) among A-type stars but somewhat contradicts 
to Dccin ct al. (2003) and Dominik & Decin (2003) who found only few young stars with 
moderate or small infrared excesses, and proposed that it might be related to the effect of 
stirring. A possible explanation of the large spread among young stars could be that the 
initial conditions of the disks (especially initial disk mass) are far from being homogeneous. 

5. Conclusions 

We searched the IRAS and ISO databases and compiled a list of debris disks exhibiting 
the highest fractional luminosity values (/^ > 10^^) in the vicinity of the Sun {d < 120 pc). 
Utilizing high-resolution far-infrared maps we attempted to exclude bogus disks from the 
sample. The fractional luminosity value for each disk was recomputed using available IRAS, 
ISO, and Spitzer data, and analysed the galactic space velocities of the objects as well as 
the distribution of the disks on the fractional luminosity versus age diagram. Our results 
are summarized as follows: 

1. We compiled a list of 60 debris disk systems of high fractional luminosity. Eleven 
of them are new discoveries and 4 out of these 11 have been confirmed by Spitzer 
observations; 

2. Disks with high fractional luminosity often belong to young stellar kinematic groups, 
providing an opportunity to obtain improved age estimates for these disks; 

3. Practically all objects with fa > 5xl0~"^ are younger than 100 Myr; 

4. The number of old systems with high f^i seems to be lower than was claimed before, 
mainly as a consequence of the age revision in connection to the young stellar kinematic 
groups; 

5. There exist many young disks of moderate fractional luminosity; 

6. Comparing the theoretical evolutionary model of Dominik & Decin (2003) with the 
observations in the fa vs. age diagram good general agreement was found. 
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A. Bogus debris disks 

In the last few years several IRAS-based debris disk candidates turned out to be bogus. 
Examples arc HD 155826 (Lisse et al. 2002), or the list of Kalas et al. (2002). Most common 
problems arc contamination by background objects (cirrus knots, galaxies), Pleiades like 
nebulosity, or unreliable point source detection by IRAS. Our list of debris disk candidates is 
based on IRAS data, however when higher resolution far-infrared observations were available 
we checked whether the above mentioned problems could have affected the detection of the 
disk. In Table 6 we fist those objects which were identified as debris system in the literature 
(and claimed to have 10~^ < fd < 10~^ in the original paper), but our analysis indicates that 
they are very likely bogus disks. In the following we briefly describe the reason of rejection. 
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HD 34739 The source position in the MIPS 70//m map differs from the star's position by 
26", but coincides with the near-infrared source 2MASS J05163646-5257397 with an offset 
of 2". 

HD 53842 This system is not a real bogus disk, since at 24//m the star shows infrared 
excess (Moor et al. 2006, in prep.) but at 70/im the IR emission comes from an independent 
compact source separated by 19". This nearby source coincides with a 2MASS J06460135- 
8359359, within a distance of 2". Due to this fact the fractional luminosity of HD 53842 
decreased below our lower limit. 

HD 56099 The source position in the MIPS 70/xm map differs from the star's position by 
24", but coincides with the near-infrared source 2MASS J07190966-F5907219 with an offset 
of 2". 

HD 72390 The peak brightness position in the ISOPHOT maps at 60 and 90/im differs 
from the stellar position by 36". Coordinates of this peak are very close to those of 2MASS 
J08 143635-8423260 which is an extended 2MASS source (XSC 1524951), with an offset of 



HD 82821 We have detected two infrared sources in the MIPS 70/im map, but their po- 
sitions significantly differ from the position of the star. The position of the brighter one 
lies at 70" from the nearby IRAS source (FSC 09319-1-0346); and is located just outside 
the 2(7 error ellipse (nearly along the major axis) but well inside the 3a error ellipse. We 
assume that this source, whose position coincides well with the near-infrared source 2MASS 
J09343630+0332417 (XSC 2391850) within a distance of 3", is responsible for the source 
confusion. The second source was probably below the IRAS sensitivity hmit. 

HD 143840 ISOPHOT mini-map observation at 90/xm and Spitzer MIPS image at lO/im 
show extended IR emission. The image of the Digitized Sky Survey also shows a reflection 
nebulosity around this star. We think that the excess far-infrared emission comes from the 
nebula. 

HD 185053 Magakian (2003) proposed that HD 185053 is the illuminating source of the 
reflection nebula GN 19.41.5. Spitzer MIPS observation at 24 and 70//m shows extended IR 
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emission around the star. We think that the excess far-infrared emission is related to the 
nebula rather than to a debris disk. 

HD 204942 The source position in the MIPS 70//m map differs from the star's position by 
22", but coincides with the near-infrared source 2MASS J21323602-2409319 with an offset 
of 4". 

HD 23484, HD 158373, HD 164330 ISOPHOT mini-map observations of these stars at 
60 and 90//m did not show excess above the photosphere. This discrepancy between the IRAS 
based excesses and the non-detection of excesses by ISOPHOT (which had better spatial 
resolution than IRAS at far-infrared wavelengths) was already mentioned by Silverstone 
(2000), who suggested cirrus contamination as the reason. 

For stars when no higher resolution data were available a set of criteria were applied to 
filter out suspicious objects which might be bogus disk (see Sect. 2.3). In Table 7 we listed 
those objects which were earlier identified as debris disk in the literature (with 10^"' < fd < 
10~^) but were rejected from the further analysis according to our criteria. Nevertheless, 
future high spatial resolution infrared data are needed to prove or disprove our judgment. 

B. New members in the young stellcir kinematic groups 

/3Pictoris moving group (BPMG). We identified two stars which are candidate mem- 
bers of this group: HD 15115 and HD 191089. HD 15115 is a northern object. Though the 
first surveys found BPMG members only in the Southern hemisphere (Barrado y Navascues 
et al. 1999; Zuckerman et al. 2001a) recently Song et al. (2003) proposed several new northern 
candidates. One of those, HIP 12545 is located within four degrees to HD 15115 on the sky. 
HD 15115 has a ROSAT counterpart with fractional X-ray luminosity of log(Lx/Lboi) = —4.94 
(logLx = 29.2 erg s~^). These properties are comparable to those of stars with similar mass 
in young star associations (see fig. 9,10 in de la Reza & Pinzon 2004). 

HD 191089 is somewhat more distant (see Table 1) than known members of the group in 
the list of Zuckerman & Song (2004b). However, large part of this list is based on a volume 
limited survey within d < 50 pc, and other authors proposed candidates at larger distances 
(Torres et al. 2002). HD 191089 shows several signs of youth. Mamajek (2004b) classified 
this star as younger than Pleiades because of its lithium abundance (EW(Li) = 95 ± 6 mA) 
was higher than 95% of Pleiades stars with similar effective temperatures. Isochronal age of 
the star (17^4 Myr, Mamajek 2004b) is in good agreement with the age of /3 Pictoris moving 
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group (12l;4Myr, see Table 4). Moreover HD 191089 rotates rapidly compared with typical 

F5 type stars {vsini = 45km s~^, Nordstrom et al. 2004). Its ROSAT X-ray luminosity 
of logLx = 29.2 erg s~^ and X-ray fractional luminosity of log(Lx/Lboi) = —4.93. These 
values are also comparable with the properties of stars with similar spectral types in nearby 
kinematic groups (see fig. 9,10 in de la Reza & Pinzon 2004). 

The GAYA2 and TucHor associations. The Great Austral Young Association2 (GAYA2) 
and Tucana/Horologium (TucHor) Association can be discussed together since they overlap 
both in their location on the projected sky plane and in velocity space. GAYA2 was dis- 
covered in the framework of the SACY (Search for Associations Containing Young stars) 
survey (Torres et al. 2002), and the known members are confined mostly in the right ascen- 
sion (R.A.) range 3^ < R.A. < 9^. Recently identified members of TucHor occupy a similar 
region (2*^ < R.A. < 7*^) and several of them show only slightly different Galactic space mo- 
tions compared to the mean UVW velocities of the Tucana nucleus and resemble the mean 
space motions of GAYA2. Although GAYA2 is more distant (located at a mean distance 
of ~ 84 pc, while TucHor members located in the same sky region have a mean distance of 
~ 50 pc), there is an overlap in radial distance, as well. Studying the relationship between 
the two associations is out of the scope of this work. As a practical solution, we assigned all 
doubtful sources (see Fig. 6) of D < 67 pc to the TucHor association and the more distant 
ones to GAYA2. Thus we propose that one of these sources, HD 37484 belongs to TucHor 
(its space velocity is not inconsistent with that of other neighbouring TucHor members). 
The measured hthium abundance of the star (Favata et al. 1993) is a strong indication of 
its youth. HD 21997, HD 30447, HD 35841, HD 38206 and HD 38207 are classified as mem- 
bers of the GAYA2 group. It is worth to mention that these five stars form a spectacular 
concentration of high fa debris disks within a relatively small area on the sky. 

Local Association. We propose that HD 10472, HD 10638, HD 218396 and HD 221853 
belong to the Local Association. HD 10472 was previously a TucHor candidate (Torres et al. 
2000; Zuckerman et al. 2001b), but recently Zuckerman &; Song (2004b) suggested that its 
membership status is uncertain, thus it may be consistent with our result. 

IC 2391 supercluster On the basis of its Galactic space velocity HD 192758 may belong 
to the 1C2391 supercluster. Its position on the CMD (see Fig. 2 and Sect. 3.2.2) also suggest 
its youth. 
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HD 110058 was earlier classified as a member of the Lower Centaurus Crux (LCC) as- 
sociation using convergent-point method (de Zeeuw et al. 1999). However, according to 
our results, its Galactic space velocity is inconsistent with the mean velocity of the LCC. 
Nevertheless HD 110058 seems to be a very young object on the basis of its position on the 
Color-Magnitude Diagram of A-type stars (see Sect. 3.2.2). 
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Table 1. Stellar properties and derived fractional luminosities of the disks 



Name Other Spectral V B-V Distance First reference Debris disk Reference fj Age Age 



Name Type to debris disk confirmation Reference 









(mag) 


(mag) 


(pc) 








(10-4) 


(Myr) 




HD 105 




GOV 


7.51 


0.595 


40 


(27) 


ISO/MIPS 


(9,27) 


2.5±0.3 




Table 4 


HD 377 




G2V 


7.59 


0.626 


40 


(9) 


MIPS 


(9) 


4.0±0.3 


[30,100] 


(4,30) 


HD 3003 




AOV 


5.07 


0.038 


46 


(29) 


MIPS* 




1.4±0.2 


•^■^-20 


Table 4 


HD 6798 




A3V 


5.60 


0.008 


83 


(20) 


IRS 


(12) 


1.6±0.3 




this work 


HD 8907 




F8 


6.66 


0.505 


34 


(27) 


ISO/MIPS 


(9,27) 


2.4±0.1 


[100,870] 


(4,35,36) 


HD 9672 


49 Get 


AlV 


5.62 


0.066 


61 


(24) 


ISO 


(33) 


9.2±0.6 


[8,20] 


(32,36) 


HD 10472 




F2IV/V 


7.62 


0.420 


67 


(27) 


MIPS* 


(21) 


3.4±0.9 


[20,150] 


Table 4 


HD 10647 




F8V 


5.52 


0.551 


17 


(31) 


ISO/MIPS* 


(6) 


3.0±0.3 


[300,7000] 


(18,36) 


HD 10638 




A3 


6.73 


0.247 


72 


(27) 






3.9±0.5 


[20,150] 


Table 4 


HD 15115 




F2 


6.79 


0.399 


45 


(27) 


ISO/MIPS 


(27) 


4.9±0.4 


12I4 


Table 4 


HD 15745 




FO 


7.47 


0.360 


64 


(27) 


ISO/MIPS 


(27) 


20.1±1.4 


<700 


(18,36) 


HD 16743 




F0/F2ni/IV 


6.78 


0.387 


60 


this work 


MIPS 




3.6±0.3 


1200l^°o 


(18) 


HD 17390 




F3IV/V 


6.48 


0.387 


45 


(27) 


ISO/MIPS* 


(1) 


1.9±0.2 


<800 


(18,36) 


HD 21997 




A3IV/V 


6.38 


0.120 


74 


(20) 


IRS 


(12) 


4.7±0.3 


^^-10 


Table 4 


HD 24966 




AOV 


6.89 


0.023 


104 


(27) 






2.4±0.5 


<100 


this work 


HD 25457 




F5V 


5.38 


0.516 


19 


(27) 


ISO/MIPS 


(9,27) 


1.0±0.2 


[50,100] 


Table 4 


HD 30447 




F3V 


7.85 


0.393 


78 


(27) 


ISO/MIPS 


(1) 


7.5±1.1 


on-l-lO 

201 10 


Table 4 


HD 32297 




AO 


8.13 


0.199 


112 


(27) 


GOR 


(26) 


33.4±6.0 


<30 


(13) 


HD 35841 




F5V 


8.91 


0.496 


104^ 


(27) 






15.5±3.7 


20l}o 


Table 4 


HD 37484 




F3V 


7.26 


0.404 


60 


(20) 


ISO/MIPS 


(9,27) 


2.7±0.5 


30I20 


Table 4 


HD 38207 




F2V 


8.47 


0.391 


IO3P 


(27) 


ISO/MIPS 


(9,27) 


10.8±0.6 




Table 4 


HD 38206 




AOV 


5.73 


-0.014 


69 


(17) 


MIPS 


(22) 


1.4±0.3 


20+"' 


Table 4 


HD 38678 


CLep 


A2Vann 


3.55 


0.104 


22 


(19) 


ISO/MIPS* 


(7) 


1.1±0.2 


200llgg 


Table 4 


HD 39060 


/3Pic 


A3V 


3.85 


0.171 


19 


(28) 


ISO/COR/SUBM 


(8,10,28) 


24.3±1.1 


12+1 


Table 4 


HD 50571 




F7ni-IV 


6.11 


0.457 


33 


this work 


MIPS* 




1.1±0.3 


ISOOlJooo 


(18) 


HD 53143 




KOIV-V 


6.81 


0.786 


18 


(17) 


ISO 


(6) 


2.0±0.5 


98011- 


(6) 


HD 54341 




AOV 


6.52 


-0.008 


93 


this work 






1.8±0.4 


<100 


this work 


HD 69830 




KOV 


5.95 


0.754 


13 


(17) 


IRS 


(3) 


2.0=* 


[600,4700] 


(29,35) 


HD 76582 




FOIV 


5.68 


0.209 


49 


(27) 






1.7±0.2 


4P;n+150 
^OU_29o 


this work 


HD 78702 




AO/AIV 


5.73 


0.000 


80 


(27) 






2.1±0.7 




this work 


HD 84870A 




A3 


7.20 


0.233 


90 


(27) 






4.5±1.0 


<520 


this work 


HD 85672 




AO 


7.59 


0.159 


93 


(27) 






4.9±1.0 


<100 


this work 


HD 92945 




KIV 


7.72 


0.873 


22 


(27) 


MIPS 


(5) 


5.3±1.2 


100 


(30,35) 



Table 1 — Continued 





Other 
NcLine 


Spsctrail 
TvDe 


V 


B-V 


Dista-iic6 
(pel 


First rsfcrsncG 
to dsbris disk 


Dsbris disk 
CO nfir iiia.t ion 


Refer ence 


(10-4) 


A ere 
(Mvrl 


A ere 
Re fere nee 


HD 107146 




G2V 


7.04 


0.604 


29 


(27) 


MIPS/COR/SUBM 


(2,9,34) 


9.2±0.9 


100+7 


(34,36) 


HD 109573A 


HR4796A 


AOV 


5.78 


0.003 


67 


(11) 


MIR/COR 


(14,25) 


47.0±2.5 




Tabic 4 


HD 110058 




AOV 


7.99 


0.148 


100 


(17) 


ISO 


(1) 


18.9±3.3 


<100 


this work 


HD 115116 




A7V 


7.07 


0.205 


85 


this work 






3.2±1.0 


<360 


this work 


HD 120534 




A5V 


7.02 


0.275 


67P 


(27) 






4.9±0.9 


<320 


this work 


HD 121812 




KO 


8.53 


0.820 


38 


this work 






14.9±4.1 


23Q+150 


this work 


HD 122106 




F8V 


6.36 


0.486 


78 


this work 






1.2±0.4 


[1000,1600] 


(18,23) 


HD 127821 




F4IV 


6.10 


0.428 


32 


(27) 


MIPS* 




2.2±0.4 


[200,3400] 


(18,36) 


HD 130693 




G6V 


8.20 


0.734 


33P 


this work 






5.8±1.9 




this work 


HD 131835 




A2IV 


7.88 


0.192 


111 


this work 






19.9±3.3 


1711 


Table 4 


HD 157728 




FOIV 


5.70 


0.229 


43 


(31) 






2.9±0.4 


<200 


this work 


HD 158352 




A8V 


5.41 


0.237 


63 


(19) 






1.8±0.5 


750ll«,° 


this work 


HD 164249 




F5V 


7.01 


0.458 


47 


(27) 


ISO/MIPS* 


(27) 


10.4±1.6 




Table 4 


HD 169666 




F5 


6.68 


0.444 


51 


this work 


MIPS 




1.5±0.3 


2200+''°° 


(18) 


HD 170773 




F5V 


6.22 


0.429 


36 


(24) 


ISO/MIPS* 


(27) 


3.8±0.4 


[200,2800] 


(18,36) 


HD 172555 




A7V 


4.78 


0.199 


29 


(19) 


MIPS* 




7.8±0.7 




Table 4 


HD 181296 


7? Tel 


AOVn 


5.03 


0.020 


48 


(17) 


ISO/MIPS* 


(27) 


2.4±0.2 


12lt 


Table 4 


HD 181327 




F5/F6V 


7.04 


0.480 


51 


(17) 


ISO/MIPS* 


(1) 


29.3±1.6 


12tl 


Tabic 4 


HD 182681 




B8/B9V 


5.66 


-0.014 


69 


this work 






1.5±0.3 


<100 


this work 


HD 191089 




F5V 


7.18 


0.480 


54 


(17) 


ISO/MIPS* 


(1) 


19.1±2.2 


12+® 


TUble4 


HD 192758 




FOV 


7.02 


0.317 


62P 


(27) 


ISO/MIPS 


(1) 


5.6±0.5 


[35,55] 


Table 4 


HD 197481 


AUMic 


MlVe 


8.81 


1.470 


10 


(27) 


MIPS/COR/SUBM 


(5,15,16) 


4.0±0.3 




Table 4 


HD 202917 




G5V 


8.65 


0.690 


46 


(27) 


MIPS* 


(27) 


2.9±0.8 




Table 4 


HD 205674 




F3/F5IV 


7.19 


0.396 


53 


this work 


MIPS* 




3.5±0.8 


2600+"°° 


(18) 


HD 206893 




F5V 


6.69 


0.439 


39 


(27) 


ISO/MIPS* 


(27) 


2.3±0.2 


<2800 


(18,36) 


HD 218396 




A5V 


5.97 


0.259 


40 


(27) 


ISO 


(27) 


2.3±0.2 


[20,150] 


Table 4 


HD 221853 




FO 


7.35 


0.405 


71 


(27) 


ISO/MIPS 


(27) 


8.0±1.1 


[20,150] 


Table 4 



Note. — Col.(l): Names. Col. (2): Other names. Cols. (3-5): Data are from the Hipparcos or the Tycho-2 Spectral Type Catalog. Col. (6): Distances. P indicates 

photometric distances otherwise Hipparcos distances were used. Col. (7): Reference for first identification as debris disks. Col. (8): Observations following the original 
IRAS discovery which independently confirmed the existence of the debris disk. COR: Coronographic observation; ISO: ISOPHOT; MIR: Observation at mid-infrared 
wavelengths; IRS: Spitzer IRS; MIPS: Spitzer MIPS; SUBM: Observation at submillimeter wavelengths. * marks those MIPS observations when we extracted only 



astromctrical information from the maps, but did not determine photometric fluxes. Col. (9): Referenees for papers related to observations of Col. (8). Col. (10): Fractional 
dust luminosities and their uncertainties. Col. (11): Ages estimates. [tl,t2] marks an age interval. Col. (12) References for age estimates. 

"■Beichman et al. (2005b) demonstrated the presence of strong spectral features for this object, therefore we took their fractional luminosity estimate rather than fitting 
the photometric points by a modified blackbody 

References. — (1) Abraham et al., 2006, in prep.; (2) Ardila ct al. (2004); (3) Bcichman ct al. (2005b); (4) Carpenter ct al. (2005); (5) Chen ct al. (2005b); (6) Decin 
et al. (2000); (7) Habing et al. (2001); (8) Heinrichscn et al.(1999); (9) http://data.spitzcr.caltcch.edu/popular/fcps/20051123-cnhanccd.vl/, FEPS Data Explanatory 
Supplement v3.0, Hines et al. (2005); (10) Holland et al. (1998); (11) Jura (1991); (12) Jura et al. (2004); (13) Kalas (2005); (14) Koerner et al. (1998); (15) Krist 
et al. (2005); (16) Liu et al. (2004); (17) Mannings & Barlow (1998); (18) Nordstrom et al. (2004); (19) Oudmaijcr et al. (1992); (20) Patten & Willson (1991); (21) 
RebuU et al. (2005); (22) Rieke et al. (2005); (23) Rocha-Pinto et al. (2004); (24) Sadakane & Nishida (1986); (25) Schneider et al. (1999); (26) Schneider et al. (2005); 
(27) Silverstone (2000); (28) Smith & Terrilc (1984); (29) Song et al. (2000); (30) Song et al. (2004); (31) Stencel & Backman (1991); (32) Thi et al.(2001); (33) Walker 
and Heinrichscn (2000); (34) WiUiams et al. (2004); (35) Wright et al. (2004); (36) Zuckerman & Song (2004a); 'Table4' refers to ages of stellar kinematic groups, as 
presented in Table 4. 
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Table 2. Summary table of fluxes measured in excess 



Name 


Instrument 


Wavelength 


Measured flux 


Flux 


Photosph. 


Ref. 








uncorrected 


corrected 


quality 


flux 








[/im] 


[Jy] 


[Jy] 




[Jy] 




HD 105 


ISOPHOT 


60 


0.107±0.014 


0.108±0.014 




3.9e-03 


1 




MIPS 


70 


0.145±0.012 


0.164±0.014 




2.9e-03 


4 




ISOPHOT 


90 


0.147±0.010 


0.159±0.011 




1.7e-03 


1 




MIPS 


160 


0.166±0.028 


0.171±0.029 




5.5e-04 


4 


HD 377 


MIPS 


24 


0.035±0.002 


0.035±0.002 




0.025 


4 




IRAS 


60 


0.179±0.014 


0.185±0.015 


3 


4.0e-03 


5 




MIPS 


70 


0.146±0.017 


0.158±0.018 




2.9e-03 


4 


HD 3003 


IRAS 


25 


0.343±0.027 


0.320±0.026 


3 


0.063 


8 


HD 6798 


IRAS 


25 


O.lllitO.018 


0.108±0.017 


3 


0.039 


8 




IRAS 


60 


0.403±0.048 


0.414±0.050 


2 


6.7e-03 


8 


HD 8907 


IRAS 


60 


0.285±0.054 


0.305±0.058 


3 


7.4e-03 


8 




ISOPHOT 


60 


0.224±0.020 


0.232±0.021 




7.4e-03 


1 




MIPS 


70 


0.232±0.007 


0.256±0.008 




5.4e-03 


4 




ISOPHOT 


90 


0.264±0.018 


0.270±0.018 




3.3C-03 


1 


HD 9672 


IRAS 


25 


0.384±0.042 


0.438±0.048 


3 


0.041 


8 




IRAS 


60 


2.02±0.12 


2.09±0.13 


3 


7.0C-03 


8 




IRAS 


100 


1.88±0.21 


1.85±0.20 


2 


2.5e-03 


8 




ISOPHOT 


120 


1.45±0.38 


1.32±0.35 




1.7e-03 


1 




ISOPHOT 


170 


0.862±0.181 


0.772±0.162 




8.7C-04 


1 


HD 10472 


IRAS 


60 


0.140±0.032 


0.145±0.033 


2 


2.4C-03 


8 


HD 10647 


IRAS 


60 


0.815±0.106 


0.877±0.114 


3 


0.024 


8 




ISOPHOT 


60 


0.803±0.040 


0.824±0.041 




0.024 


1 




IRAS 


100 


1.08±0.17 


1.08±0.17 


2 


8.7e-03 


8 


HD 10638 


IRAS 


60 


0.348±0.038 


0.360±0.040 


3 


3.6e-03 


8 


HD 15115 


MIPS 


24 


0.055±0.006 


0.054±0.006 




0.032 


7 




IRAS 


60 


0.441±0.048 


0.473±0.052 


3 


5.1e-03 


8 




ISOPHOT 


60 


0.401±0.019 


0.415±0.020 




5.1e-03 


1 




ISOPHOT 


90 


0.419±0.029 


0.427±0.030 




2.2e-03 


1 


HD 15745 


MIPS 


24 


0.162±0.016 


0.168±0.017 




0.016 


7 




IRAS 


25 


0.177±0.028 


0.207±0.033 


2 


0.015 


8 




IRAS 


60 


0.868±0.061 


0.875±0.061 


3 


2.6e-03 


8 




ISO 


60 


0.753±0.038 


0.770±0.039 




2.6e-03 


1 




ISO 


90 


0.532±0.037 


0.515±0.036 




l.lc-03 


1 


HD 16713 


MIPS 


24 


0.019±0.005 


0.019±0.005 




0.029 


7 




IRAS 


()0 


().:i()l±0.():i:-! 


0.:-!2:-!±0.0:-!() 


O 


l.(ic-0:-! 


8 


HD 17390 


IRAS 


60 


0.226±0.034 


0.240±0.036 


3 


6.4e-03 


8 




ISO 


60 


0.253±0.021 


0.262±0.022 




6.4e-03 


1 




ISO 


90 


0.228±0.023 


0.231±0.023 




2.8e-03 


1 


HD 21997 


IRAS 


60 


0.595±0.036 


0.646±0.039 


3 


3.9e-03 


8 




IRAS 


100 


0.636±0.108 


0.629±0.107 


2 


1.4e-03 


8 


HD 24966 


IRAS 


60 


0.184±0.031 


0.190±0.032 


3 


1.9e-03 


8 


HD 25457 


MIPS 


24 


0.202±0.011 


0.204±0.011 




0.147 


4 




IRAS 


60 


0.287±0.075 


0.285±0.074 


2 


0.023 


8 




ISOPHOT 


60 


0.293±0.030 


0.297±0.030 




0.023 


1 
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Table 2 — Continued 



Name 


Instrument 


Wavelength 


Measure 


d flux 


Flux 


Photosph. 


Ref. 








uncorrected 


corrected 


quality 


flux 








[/urn] 


[Jy] 


[Jy] 




[Jy] 






MIPS 


70 


0.288±0.012 


0.310±0.013 




0.017 


4 




ISOPHOT 


90 


0.246±0.018 


0.242±0.018 




0.010 


1 


HD 30447 


MIPS 


24 


0.028±0.003 


0.028±0.003 




0.012 


7 




IRAS 


60 


0.274±0.030 


0.296±0.033 


3 


1.9C-03 


8 




ISO 


60 


0.247±0.062 


0.256±0.064 




1.9e-03 


1 




ISO 


90 


0.271±0.065 


0.277±0.067 




8.3e-04 


1 


HD 32297 


IRAS 


25 


0.211±0.034 


0.256±0.041 


2 


5.7e-03 


8 




IRAS 


60 


1.12±0.07 


1.14±0.07 


3 


9.7e-04 


8 


HD 35841 


IRAS 


60 


0.188±0.040 


0.195±0.041 




8.1C-04 


8 


HD 37484 


MIPS 


24 


0.053±0.003 


0.055±0.003 




0.021 


4 




IRAS 


60 


0.128±0.036 


0.124±0.035 


2 


3.3e-03 


8 




ISOPHOT 


60 


0.094±0.012 


0.095±0.012 




3.3e-03 


1 




MIPS 


70 


0.104±0.011 


0.110±0.012 




2.4e-03 


4 




ISOPHOT 


90 


0.047±0.012 


0.044±0.011 




1.5e-03 


1 


HD 38207 


MIPS 


24 


0.017±0.003 


0.017±0.003 




6.9e-03 


4 




ISOPHOT 


60 


0.205±0.015 


0.212±0.015 




l.le-03 


1 




MIPS 


70 


0.178±0.006 


0.194±0.007 




7.9e-04 


4 




ISOPHOT 


90 


0.177±0.013 


0.177±0.013 




4.8e-04 


1 


HD 38206 


MIPS 


24 


0.115±0.023 


0.119±0.024 




0.033 


9 




IRAS 


25 


0.110±0.021 


0.113±0.021 


2 


0.030 


8 




IRAS 


60 


0.313±0.038 


0.313±0.038 


3 


5.1e-03 


8 


HD 38678 


IRAS 


25 


1.16±0.09 


0.97±0.08 


3 


0.296 


8 




IRAS 


60 


0.366±0.073 


0.302±0.060 


3 


0.051 


8 




ISOPHOT 


60 


0.404±0.020 


0.390±0.019 




0.051 


1 


HD 39060 


IRAS 


12 


3.40±0.14 


2.38±0.10 


3 


1.09 


8 




IRAS 


25 


8.81±0.35 


10.29±0.41 


3 


0.254 


8 




IRAS 


60 


19.7±0.79 


18.92±0.76 


3 


0.043 


8 




IRAS 


100 


ll.U±0.44 


y.32±0.37 


3 


0.015 


8 




ISOPHOT 


120 


7.9itl.2 


6.68±1.01 




0.011 


1 




ISOPHOT 


150 


4.80±0.73 


4.39±0.67 




6.8e-03 


1 




ISOPHOT 


170 


4.0±0.6 


3.29±0.49 




5.4e-03 


1 




ISOPHOT 


200 


2.10±0.35 


2.03±0.34 




3.9e-03 


1 


HD 50571 


IRAS 


60 


0.183±0.040 


0.186±0.041 


2 


0.011 


8 


HD 53143 


IRAS 


60 


0.152±0.040 


0.155±0.040 


2 


0.012 


8 




ISOPHOT 


60 


0.214±0.011 


0.219d=0.011 




0.012 


1 




ISOPHOT 


90 


0.145±0.011 


0.142±0.011 




5.1e-03 


1 


HD 54341 


IRAS 


60 


0.203±0.047 


0.210±0.048 


2 


2.7e-03 


8 


HD 69830 


IRAS 


25 


0.341±0.038 


0.292±0.032 


3 


0.142 


8 


HD 76582 


IRAS 


60 


0.391±0.047 


0.403±0.048 


3 


9.2e-03 


8 


HD 78702 


IRAS 


60 


0.314±0.047 


0.327±0.049 


2 


5.6C-03 


8 




IRAS 


100 


0.954±0.267 


0.987±0.276 


2 


2.0e-03 


8 


HD 84870A 


IRAS 


60 


0.243±0.046 


0.252±0.048 


3 


2.4e-03 


8 


HD 85672 


IRAS 


60 


0.193±0.037 


0.200±0.038 


3 


1.4C-03 


8 


HD 92945 


IRAS 


60 


0.248±0.045 


0.270±0.049 


3 


6.3e-03 


8 
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Table 2 — Continued 



Name 


Instrument 


Wavelength 


Measured flux 


Flux 


Photosph. 


Ref. 








uncorrected 


corrected 


quality 


flux 








[fj,m] 


[Jy] 


[Jy] 




[Jy] 






Riirs 


70 


().271±0. ()()() 


l).:-!01±0.l)()7 




l.lie-O:-! 


2 


HD 107146 


MIPS 


24 


0.060±0.003 


0.059±0.003 




0.042 


4 




IRAS 


60 


0.705±0.056 


0.777±0.062 


3 


6.7e-03 


8 




MIPS 


70 


0.648±0.014 


0.727±0.016 




4.9e-03 


4 




IRAS 


100 


0.910±0.155 


0.913±0.155 


2 


2.4e-03 


8 




SCUBA 


450 


0.130±0.040 


0.130±0.040 




l.le-04 


10 




SCUBA 


850 


0.020±0.004 


0.020±0.004 




3.2e-05 


10 


HD 109573A 


IRAS 


12 


0.284±0.011 


0.214±0.009 


3 


0.131 


5 




IRAS 


25 


3.73ib0.00^ 


4.38ib0.00^ 


3 


0.031 


5 




IRAS 


60 


o.U 1 ±U.UU 


i . ( IzLU.UU 


3 


5.2e-03 


5 




IRAS 


100 


4.U9d=U.U8 


3.8o±U.U8 


3 


1.9e-03 


5 




SCUBA 


850 


n ni Q-Un r\f\A 
U.UiyitU.UU4 


n ni Q-un nn/i 




2.5e-05 


3 


TTT\ "1 1 f\f\fCi 

HD 110058 


IRAS 


25 


0.2oD±0.029 


0.300±0.033 


3 


5.7e-03 


8 




IRAS 


60 


U.oOCSlllU.UDZ 


n ^41 -i-H nf^s 


2 


9.8e-04 


8 




io^Jr^xltJ ± 


OU 


u.oo7±U.Uoo 


u.oo7±U.Uoo 




y.oe-U4 


i 




IbUFHL) 1 


90 


U.^ioitU.UoU 


n oni -Un n/i fi 




4.3e-04 


1 


HD 115116 


IRAS 


60 


n 1 noj_n c\f^A 
U.iyziLU.Uo4 


n 1 riQ-un nc^c 


2 


2.5e-03 


8 


TTT\ "1 fir^^O A 

HD 120534 


IRAS 


60 


f\ IIP' 1 r\ f\ A ^ 


0.32d±0.049 




3.3e-03 


8 


HD 121812 


IRAS 


60 


U.2do±U.U47 


U.2oy±U.Uo2 


3 


2.7e-03 


8 




IHAb 


100 


f\ A a A 1 

0.484±0.135 


0.492±0.138 


2 


9.6e-04 


8 


HD 122106 


IRAS 


60 


0.161±0.042 


0.164±0.043 


2 


9.2e-03 


8 


HD 127821 


IRAS 


60 


O.o44it0.041 


0.37d±0.045 


3 


0.011 


8 




IRAS 


100 


0.525±0.137 


0.529±0.138 


2 


3.8e-03 


8 


HD 130693 


IRAS 


60 


0.149±0.042 


0.154±0.043 




3.2e-03 


8 


HD 131835 


IRAS 


25 


0.186±0.034 


0.224±0.040 


2 


6.1e-03 


8 




IRAS 


60 


0.684±0.062 


0.681±0.061 


3 


l.Oe-03 


8 


HD 157728 


IRAS 


25 


0.213±0.019 


0.221±0.020 


3 


0.053 


8 




IRAS 


60 


0.536±0.048 


0.531±0.048 


3 


9.0e-03 


8 


HD 158352 


IRAS 


25 


0.258±0.044 


0.255±0.043 


3 


0.075 


8 




IRAS 


60 


0.366±0.062 


0.348±0.059 


2 


0.013 


8 


HD 164249 


IRAS 


60 


0.647±0.091 


0.669±0.094 


2 


4.7e-03 


8 




ISOPHOT 


60 


0.740±0.037 


0.761±0.038 




4.7e-03 


1 




ISOPHOT 


90 


0.568±0.040 


0.560±0.039 




2.1e-03 


1 


HD 169666 


MIPS 


24 


0.087±0.009 


0.090±0.009 




0.039 


7 




IRAS 


2.1 


O.lDliO.Dl:-! 


0.()88±0.()11 




o.o:«) 


S 


HD 170773 


IRAS 


60 


0.504±0.086 


0.555±0.094 


2 


9.1e-03 


8 




ISOPHOT 


60 


0.547±0.027 


0.562±0.028 




9.1e-03 


1 




ISOPHOT 


90 


0.771±0.054 


0.821±0.057 




4.0e-03 


1 


HD 172555 


IRAS 


12 


1.52±0.06 


1.34±0.05 


3 


0.508 


8 




IRAS 


25 


1.09±0.06 


0.91±0.05 


3 


0.119 


8 




IRAS 


60 


0.306±0.046 


0.241±0.036 


3 


0.020 


8 


HD 181296 


IRAS 


12 


0.481±0.024 


0.351±0.018 


3 


0.264 


8 




IRAS 


25 


0.491±0.025 


0.496±0.025 


3 


0.061 


8 




IRAS 


60 


0.495±0.045 


0.449±0.040 


3 


0.011 


8 
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Table 2 — Continued 



Name 


Instrument 


Wavelength 


Measured flux 


Flux 


Photosph. 


Ref. 








uncorrected 


corrected 


quality 


flux 








[fj,m\ 


r T, 1 

|JyJ 


|JyJ 




[JyJ 






TSOPHOT 


60 


0.436±0.022 


0.433±0.022 




0.011 


1 




TSOPHOT 


90 


0.307±0.022 


0.286±0.021 




4.6c-03 


1 


HD 181327 


IRAS 


25 


0.248±0.020 


0.286±0.023 


3 


0.027 


8 




IRAS 


60 


1.86±0.07 


1.93±0.08 


3 


4.7e-03 


8 




ISOPHOT 


60 


1.69±0.17 


1.73±0.17 




4.7e-03 


1 




ISOPHOT 


90 


1.40±0.14 


1.41±0.14 




2.1e-03 


1 




IRAS 


100 


1.72±0.21 


1.69±0.20 


2 


1.7e-03 


8 




ISOPHOT 


170 


0.820±0.214 


0.736±0.192 




5.9e-04 


1 


HD 182681 


IRAS 


60 


0.412±0.066 


0.426±0.068 


2 


5.6C-03 


8 


HD 191089 


IRAS 


25 


0.339±0.051 


0.387±0.058 


2 


0.024 


8 




IRAS 


60 


0.711±0.057 


0.729±0.058 


3 


4.1e-03 


8 




ISOPHOT 


60 


0.774±0.080 


0.781±0.081 




4.1e-03 


1 




ISOPHOT 


90 


0.394±0.040 


0.370±0.038 




1.8e-03 


1 


HD 192758 


MIPS 


24 


0.042±0.004 


0.041±0.004 




0.022 


7 




IRAS 


60 


0.360±0.050 


0.388±0.054 




3.5e-03 


8 




ISOPHOT 


60 


0.370±0.022 


0.383±0.023 




3.5e-03 


1 




ISOPHOT 


90 


0.375±0.026 


0.384±0.027 




1.6e-03 


1 


HD 197481 


IRAS 


60 


fi 269+0 046 


yj. ^tJ^_l_\J.\J'-±tJ 


2 


0.024 


8 




MIPS 


70 


0.196±0.015 


0.207±0.016 




0.017 


2 




SHARC 


350 


0.072±0.020 


0.072±0.020 




6.8e-04 


2 




SCUBA 


850 


0.014±0.002 


0.014±0.002 




l.le-04 


6 


HD 202917 


ISOPHOT 


60 


0.046±0.014 


0.047±0.014 




1.9e-03 


1 


HD 205674 


IRAS 


60 


0.199±0.040 


0.205±0.041 


3 


3.7e-03 


8 


HD 206893 


IRAS 


60 


0.228±0.039 


0.247±0.042 


3 


6.3C-03 


8 




ISOPHOT 


60 


0.234±0.015 


0.242±0.016 




6.3C-03 


1 




ISOPHOT 


90 


0.259±0.018 


0.268±0.019 




2.8e-03 


1 


HD 218396 


IRAS 


60 


0.410±0.061 


0.450±0.068 


3 


8.6e-03 


8 




ISOPHOT 


60 


0.400±0.020 


0.412±0.021 




8.6e-03 


1 




ISOPHOT 


90 


0.553±0.039 


0.585±0.041 




3.8e-03 


1 


HD 221853 


MIPS 


24 


0.079±0.008 


0.082±0.008 




0.019 


7 




IRAS 


60 


0.327±0.062 


0.329±0.062 


3 


3.0C-03 


8 




ISO 


60 


0.368±0.031 


0.376±0.032 




3.0e-03 


1 




ISO 


90 


0.231±0.016 


0.223±0.015 




1.3e-03 


1 



Note. — Col.(l): Names. Col. (2): Instrument. Col. (3): Wavelength. Cols. (4-5): Measured flux density and uncertainty at 
the specific wavelength. On the infrared flux density values in Col. (5) color correction was applied. No color correction was 
applied on submillimeter fluxes. CoI.(6): Quality of flux density, if available. Col.(7): Predicted flux density of the stellar 
photosphere at the specific wavelength. For a detailed description of the prediction method, see Sect. 2.1. Col. (8): References 
for papers related to the quoted fiux and its uncertainty in Cols. (4-5). 

^No quoted fiux uncertainty in the IRAS Serendipitous Survey Catalog. We assumed SF25 = 0.15 Jy, and SFgo = 0.32 Jy 
(4% relative uncertainty in both cases). 

References. — (1) Abraham et al., 2006, in prep.; (2) Chen et al. (2005b); (3) Greaves et al. (2000); (4) 
http://data.spitzer.caltech.edu/populaT/feps/20051123-enhanced-vl/, FEPS Data Explanatory Supplement v3.0, Hines et al. 
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(2005); (5) Klcinmann et al. (1986); (6) Liu et al. (2004); (7) Moor et al. (2006), in prep.; (8) Moshir et al. (1989); (9) Rieke 
et al. (2005); (10) Williams et al. (2004). 
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Table 3. Radial velocity information and membership in moving groups 



Name 


Vrad 

n —11 
[km s 


Ref. 


U 

n —11 
[km s ^\ 


V 

n —11 
[km s ^] 


W 

n —11 
[km s ^] 


Assoc. SKG 


Prob. 

[%] 


Ref. 


HD 


105 


+1.6 ± 0.3 


7 


—9.8 ± 0.4 


— 21.5 ± 0.8 


— 1.3 ± 0.3 


TiirHor 


78 


Q 


HD 


377 


+1.3 ± 0.3 


7 


— 14.4 + 0.6 


—7.0 ± 0.4 


—3.9 ± 0.3 






15 


HD 


3003 


+7.0 ± 2.0 


13 


—9.6 + 0.8 


—20.9 ± 1.0 


—0.6 ± 1.6 




85 


15 


HD 


6798 


+ 10 + 43 


5 


—36 6 + 2 7 


— 13 5 + 3 5 


+4.1 ± 1.3 








HD 


8907 


+8.8 ± 0.4 


7 


—10.4 ± 0.3 


—4.6 ± 0.4 


— 16.7 ± 0.4 








HD 


9672 


+11.4 ± 1.8 


5 


—23.9 ± 1.1 


— 16 6 + 08 


—6.6 ± 1.7 








HD 


10472 


11 Q C -L 1 K 

-L C ij _l_ -L.*J 




—8.5 ± 0.7 


—24.2 + 1.1 


— 10.4 ± 1.3 


LA 


38 


UlllD W U± iV 


HD 


10647 




7 


— 1.2 ± 0.1 


_2R fi + 2 


— 17.8 ± 0.2 








HD 


10638 


—0.4 ± 1.2 


11 


— 12.3 ± 1.0 


— 25 + 1 5 


— 14.8 + 1.0 


LA 


29 


this work 


HD 


15115 




Llilo WWiJV 


— 13.2 ± 1.9 


— 17.8 ± 1.2 


—6.0 ± 2.3 


BPMG 


22 


Ulllo W \JL Iv 


HD 


15745 


+ 10 5 + 12 


■J-Viic worlc 


— 16.5 ± 1.1 


— 10.8 ± 1.3 


— 10.7 ± 0.7 








HD 


16743 


+21.9 ± 1.1 


this work 


—23 6 + 09 


— 18 + 06 


— 15.1 ± 0.9 








HD 


17390 


+7.2 ±1.8 


7 


_i K Q _i_ n 8 

i-Kj .-J —1— \J . O 


—9.6 + 0.5 


+1.0 ±1.6 








HD 


21997 


+ 17.3 ± 0.8 


5 


— 12 9 + 05 


— 22 3 + 1 


—3.9 ± 0.9 


GAYA2 


62 


tliis work 


HD 


24966 


+29 6 + 2 2 


til is work 


— 15 + 09 


— 26.4 + 1.4 


— 13.2 ± 1.8 








HD 


25457 


4_i 7 « -1- n 2 


7 


—7.9 ± 0.2 


—28.7 ± 0.4 


— 11.9 ± 0.2 


ABDor 


30 


16 


HD 


30447 


+21 3 + 25 


vixio vv yj\- IV 


—13.1 ± 1.4 


—20 9 + 1 6 


—3.8 ± 1.7 


GAYA2 


54 


uxiio w \j 1. IV 


HD 


32297 


1 21 Q + 2 1 


"1" ni G TXTTti'l;' 


_if; Q -1- 2 

XiJ.tJ _1_ ^.\J 


—16.0 ± 1.5 


—11.0 ± 1.0 








HD 


35841 


+23 1 + 13 


this work 


—13.0 ± 1.0 


—21.3 ± 1.7 


—3.6 ± 1.7 


GAYA2 


56 


this work 


HD 


37484 


+23 + 2 6 


this work 


— 11.6 ± 1.4 


—20.4 ± 1.9 


—5.2 ± 1.3 


TucHor 


25 


this work 


HD 


38207 


+24.9 ± 1.4 


LlllO W \JV xv 


— 14.7 ± 1.1 


—21.2 ± 1.5 


—4.1 ± 1.5 


GAYA2 


32 




HD 


38206 


+24.9 ± 0.6 


5 


— 13.7 ± 0.5 


—21.2 ± 0.5 


—6.1 + 0.4 


GAYA2 


27 


this work 


HD 


38678 


11 o Q -1- 2 7 


5 


_i Q fi -1- 2 


— 10.5 ± 1.6 


—8.1 ± 1.0 


v_'cio \j\jy 


28 


\ 


HD 


39060 


+20 2 + 04 


14 


—10.9 ± 0.1 


—16 2 + 03 


—9.2 ± 0.2 


BPMG 


98 


2 


HD 


50571 


+22.2 ± 1.2 


7 


—16.2 ± 0.4 


—22.3 ± 1.1 


—4.4 ± 0.5 








HD 


53143 


+21 9 + 01 


7 


—25 5 + 03 


— 18.1 ± 0.1 


— 15.2 ± 0.1 








HD 


54341 


+47.4 + 1.4 


this work 


— 16 5 + 05 


—44.0 ± 1.3 


—8.7 ± 0.5 








HD 


69830 


1 2Q 8 + 01 


7 


128 Q + 5 


—60 9 + 04 


— 10.1 ± 0.2 








HD 


76582 


—4.0 ±2.5 


11 


+10.7 ±1.8 


+5.1 ± 1.1 


+9.9 ± 1.5 








HD 


78702 


+16.9 ± 2.1 


5 


-27.5 ± 1.6 


-10.5 + 1.8 


-7.3+1.1 








HD 


84870A 


+3.7 + 0.4 


this work 


-7.8 ±0.5 


-25.5 + 2.0 


-5.7 ±0.7 








HD 


85672 


-5.2 ±6.8 


4 


-9.5 ±4.2 


-3.0+1.6 


-14.8 ±5.4 








HD 


92945 


+22.6 + 0.2 


7 


-15.2 ±0.3 


-27.8 + 0.2 


-4.3 ±0.3 








HD 


107146 


+1.5 + 0.2 


7 


-10.6 ±0.3 


-28.8 + 0.7 


-5.2 + 0.3 








HD 


109573A 


+9.4 ±2.3 


5 


-9.0 ± 1.3 


-19.0 + 1.9 


-4.4+1.0 


TWA 


64 


10 


HD 


110058 


+21.7+1.3 


this work 


+0.0 + 1.3 


-26.6 ± 1.4 


-2.1 + 0.8 








HD 


115116 


-2.3 + 0.7 


4 


-30.5 + 1.5 


-17.2 + 1.1 


-1.8 + 0.3 








HD 


120534 


+46.7+1.1 


this work 


+21.8 + 1.8 


-43.1 ±3.2 


+15.0+1.8 








HD 


121812 


-15.2 ±0.7 


8 


-29.6 + 1.2 


-61.3 ±2.7 


-1.8 + 0.9 








HD 


122106 


-1.6 ±5.0 


7 


+0.6 + 2.6 


-21.7 ±2.0 


-9.0 + 4.2 








HD 


127821 


-15.4 ±2.6 


7 


-16.0 + 0.5 


-26.2 ± 1.6 


-2.7 + 2.0 








HD 


130693 


+14.3 + 0.8 


this work 


+8.1 + 0.9 


-7.3 ±0.6 


+10.4 + 0.8 








HD 


131835 


+3.3 ± 1.7 


this work 


-4.6 + 1.6 


-17.2 ±2.1 


-3.0 + 1.0 


UCL 


50 


12 


HD 


157728 


-19.7 ± 1.2 


11 


-7.0 ± 0.8 


-21.3 ±0.8 


-4.4 ±0.6 








HD 


158352 


-36.1 + 1.2 


11 


-35.5 + 1.1 


-19.9 ±0.5 


+7.4 + 0.9 
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Table 3 — Continued 



Name 


V.ad 


Ref. 


U 


V 


W 


Assoc. SKG 


Prob. 


Ref. 




[km s-l] 




[km s-l] 


[km s-i] 


[km s-i] 




[%] 




HD 164249 


-0.2 ±0.5 


7 


-7.6 + 0.6 


-15.3 + 0.7 


-8.9 + 0.4 


BPMG 


20 


9 


HD 169666 


-44.3 ±0.6 


7 


+1.3 + 0.3 


-46.8 + 0.6 


-7.9 + 0.5 








HD 170773 


-26.3 ± 1.1 


7 


-30.5 + 1.1 


-3.8 + 0.2 


-12.3 ± 0.6 








HD 172555 


+2.0 ±2.5 


11 


-11.0 + 2.0 


-15.6 + 1.2 


-9.3 + 1.0 


BPMG 


96 


9 


HD 181296 


-2.0 ± 10.0 


13 


-10.7 + 8.6 


-14.9 + 2.7 


-7.3 + 4.4 


BPMG 


87 


9 


HD 181327 


+0.2 ±0.4 


7 


-9.1 + 0.5 


-16.2 + 0.7 


-8.4 + 0.4 


BPMG 


57 


9 


HD 182681 


+1.4 + 5.0 


3 


-0.4 + 4.6 


-13.4+1.1 


-10.8+1.8 








HD 191089 


-5.8 ± 2.2 


7 


-7.8 + 1.9 


-16.2 + 0.9 


-10.3 ± 1.2 


BPMG 


38 


this work 


HD 192758 


-11.1 ± 1.2 


this work 


-18.4 + 2.1 


-13.8 + 2.8 


-6.7 + 2.7 


IC2391 


50 


this work 


HD 197481 


-4.5 ± 1.3 


5 


-10.1 + 1.0 


-16.4 + 0.3 


-10.5 + 0.8 


BPMG 


56 


2 


HD 202917 


-1.6 ±0.2 


7 


-8.2 + 0.4 


-20.0+1.1 


-0.3 + 0.2 


TucHor 


48 


15 


HD 205674 


+1.1 + 5.1 


this work 


-1.6 + 3.0 


-24.5 + 2.4 


-16.6 + 3.7 








HD 206893 


-12.9 + 1.4 


7 


-20.0 + 0.9 


-7.7 + 0.7 


-2.1 + 1.0 








HD 218396 


-12.6 + 1.3 


5 


-12.4 + 0.5 


-21.4 + 1.1 


-7.4 + 0.9 


LA 


62 


this work 


HD 221853 


-4.2 + 2.1 


this work 


-12.7 + 0.9 


-20.6 + 1.8 


-11.2 + 1.9 


LA 


94 


this work 



Note. — Col.(l): Names. Col. (2): Radial velocity and its uncertainty. Col. (3): Reference for the source of measurement. 
Cols. (4-6): U,V,W Galactic space velocity components, U is positive towards the Galactic centre, V is positive in the direction 
of galactic rotation and W is positive towards the North galactic pole. Col. (7): Assigned stellar kinematic group (see Table4). 
Col. (8): Membership probability. Col. (9); Reference for membership identification. 

References. — (1) Barrado y Navascucs et al. (1998); (2) Barrado y Navascucs et al. (1999); (3) Evans (1967); (4) Grenier et 
al. (1999); (5) Kharchenko et al. and references therein (2004); (6) Mamajek et al. (2004a); (7) Nordstrom et al. and references 
therein (2004); (8) Strassmeier et al. (2000); (9) Song et al. (2003); (10) Webb et al.(1999); (11) Wilson (1953); (12) de Zeeuw 
et al. (1999); (13) Zuckerman & Webb (2000); (14) Zuckerman et al. (2001a); (15) Zuckerman et al. (2001b); (16) Zuckerman 
& Song (2004b) 
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Table 4. Description of stellar kinematic groups 



Group Name 


U,V,W 


Ref. 


Age 


Ref. 




[km s-i] 




[Myr] 




AB Dor Moving Group 


-8, -27, -14 


13 


50 


13 








100 


3 


/3Pictoris Moving Group 


-11, -16, -9 


13 


12+_1 


11 


Castor Moving Group 


-10.7, -8.0, -9.7 


1 


200 ± 100 


1 


Great Austral Young Association2 


-11.0, -22.5, -4.6 


9 


20 


9 


Hyades cluster 


-40, -17, -3 


13 


600 


13 


IC 2391 supercluster 


-20.6, -15.7, -9.1 


5 


35-55 


5 


Local Association 


-11.6, -21.0, -11.4 


5 


20-150 


5 


Lower Centaurus Crux 


-8.2, -18.6, -6.4 


6 


16 ± 1 


4 


Upper Centaurus Lupus 


-6.8, -19.3, -5.7 


6 


17 ± 1 


4 


Tucana /Horologium 


-11, -21, 


13 


30 


13 








10-30 


7 








20 


8 








10-40 


12 


TW Hydrac Association 


-11, -18, -5 


13 


8 


13 








5-15 


10 


Ursa Major Moving Group 


4-14, -1-1, -9 


13 


300 


13 








500 


2 



References. — (1) Barrado y Navascucs ct al. (1998); (2) King et al. (2003); (3) 
Luhman et al. (2005); (4) Mamajek et al. (2002); (5) Montes et al. and references 
therein (2001); (6) Sartori et al. (2003); (7) Stelzer & Neuhauser (2000); (8) Torres 

ct al. (2000): (9) Torres ct al. (2002); (10) Weintraub ct al. (2000): (11) Zuckcrman 
et al. (2001a); (12) Zuckerman et al. (2001b); (13) Zuckcrman & Song and references 
therein (2004b); 
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Table 5. Comparison of age estimates for newly identified moving group members 



Name 


Age estimates 
[Myr] 


References 


HD 10472 


2000+1°°° 
30 
[20,150] 


5 
8 

this study 


HD 15115 


900+ii^°° 
100 


5 
8 

this study 


HD 21997 


100 
20+1° 


8 

this study 


HD 30447 


2100+™° 
< 100 
20+1° 


5 
8 

this study 


HD 37484 


1500+«°° 

— ioUU 

[30,100] 
30+1° 


5 
1 

this study 


HD 38207 


[100,300] 

20+1° 


1 

this study 


HD 38206 


~ 9+1^ 

2oi;° 


3 

this study 


HD 191089 


17±« 
30001^°° 

< 100 

12+_l 


4 
5 

8 

this study 


HD 218396 


'■J^-682 

30 
[20,150] 


7 
8 

this study 


HD 221853 


800 
1800 
1700+5°° 
< 100 
[20,150] 


6 
2 
5 
8 

this study 



References. — (1) Carpenter ct al. (2005); 
(2) Decin et al. (2000); (3) Gerbaldi et al. 
(1999); (4) Mamajek (2004); (5) Nordstrom et 

al. (2004); (6) Silvcrstonc (2000); (7) Song et 
al. (2001); (8) Zuckerman & Song (2004a); 
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Table 6. List of bogus disks 



Name 


Nearby IRAS 


First reference 


Reason 


Instrument 


Position of unrelated 






Source 


as debris disk 


of rejection 




nearby IR source 








candidate 






RA(2000) 


DEC (2000) 


HD 


23484 


F03423-3826 


2& 


no detectable excess'' 


ISOPHOT 






HD 


34739 


F05154-5301 


2 


source confusion 


MIPS 


5 16 36.3 


-52 57 40 


HD 


53842 


F06539-8355 


1 


source confusion'' 


MIPS 


6 46 02.8 


-83 59 37 


HD 


56099 


F07149+5913 


2 


source confusion 


MIPS 


7 19 09.6 


-1-59 07 20 


HD 


72390 


F08210-8414 


2 


source confusion 


ISOPHOT 


08 14 39 


-84 23 23 


HD 


82821 


F09319+0346 


2 


source confusion 


MIPS 


09 34 36.2 


+03 32 39 


HD 


143840 


F16001-0440 


2 


extended emission 


ISOPHOT, MIPS 






HD 


158373 


F17265-0957 


23. 


no detectable excess" 


ISOPHOT 






HD 


164330 


F17559+6236 


2& 


no detectable excess" 


ISOPHOT 






HD 


185053 


F19415-8123 


2 


extended emission 


MIPS 






HD 


204942 


F21297-2422 


2 


source confusion 


MIPS 


21 32 35.8 


-24 09 30 



Note. — Col.(l): Name. Col. (2): Identification of nearby IRAS source. The IRAS source is always located within 30" of 
the star position. Col. (3); Reference for first mention as debris disk. Col. (4): Reason why the object was classified as bogus 
disk and rejected from further analyses. Col. (5); Instrument. Cols. (6- 7): When the object was nominated as bogus disk due 
to source confusion, position of the unrelated nearby IR source is given. In the cases of MIPS observations MOPEX was 
used to extract source coordinates from 70/im MIPS images. The coordinates of peak brightness in ISOPHOT maps were 
determined by fitting a point source jirofile. 

"Silverstone (2000) selected as debris disk candidate on the basis of IRAS data, but found that ISOPHOT observations 
did not confirm the detection of IR excess. He proposed that the non-detection of far-infrared flux excess towards the star's 
position can be explained by cirrus contamination. 

''HD 53842 is not really bogus. At 24/im it shows IR excess above the photosphere (see Moor et al. 2006, in prep.). 
However at 70//m the excess emission is related to a nearby infrared source. 



References. — (1) Mannings & Barlow (1998); (2) Silverstone (2000). 
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Table 7. List of rejected suspicious objects 



Name 


Nearby IRAS 


First reference 


Reason of 


Name of 




Source 


as debris disk 
candidate 


suspicion 


background source 


HIP 13005 


F02444+1505 


2 


nearby extended 2MASS source'' 


2MASS J02471368+1518315 (XSC 122165) 


HD 33095 


F05049-1927 


2 


cirrus 




HD 36162 


F05275+1519 


2 


nearby extended 2MASS source 


2MASS J05303093+1521513 (XSC 2620748) 


HD 39944 


F05526-2535 


1 


nearby galaxy'' 


ESO 488-41 


HD 83870 


F09393+4111 


2 


nearby galaxy 


PGC 27759 


HD 97455 


F11107+5541 


2 


nearby galaxy 


PGC 34197 


HD 124718 


F14129-2707 


2 


nearby 2MASS source with an 
excess in the Kg band 


2MASS J14155109-2721050 


HD 140775 


F15429+0536 


2 


star locates in tlie wall of 
Local Bubble or beyond 




HD 154145 


F17011-0004 


2 


star locates in the wall of 
Local Bubble or beyond 





Note. — Col.(l); Name. Col.(2): Identification of nearby IRAS source. The IRAS source is always located within 30" of the star 
position. CoL(3): Reference for first mention as debris disks. Col. (4): Reason of suspicion. For more details see Sect. 2.3. Col. (5): 

Name of background source. 

^Zuckerman & Song (2004a) also found suspicious this debris candidate based on the offset between its IRAS positions measured 
at 12 and 60/um. 

''Sylvester & Mannings (2000) also noted this coincidence between the position of the IRAS source and the nearby galaxy. 



References. — (1) Mannings & Barlow (1998); (2) Silverstone (2000). 
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Fig. 1. — U,V and W,V planes for stars in Table 3. The dashed rectangle marks the young 
disk population as defined by Leggett (1992). 
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Fig. 2. — Color-magnitude diagram of A-typc stars. Objects measured by Hipparcos within 
lOOpc from the Sun with parallax error (<10%) and (B-V) error (<0.01 mag) are represented 
by small dots. Filled symbols mark positions of A- and FO-type stars from Table 1; triangles 
correspond to objects assigned to any stellar kinematic group in this work. Stars of a Per 
and Hyades open clusters are represented by squares and asterisks, respectively. 
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Fig. 3. — Histogram of distances in the 3D velocity space to the closest moving group. Solid 
line: stars from Table 3; dashed line: the volume-limited Hipparcos sample (see Sect. 4.1). 
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Fig. 4. — Fractional luminosity of the infrared excess as a function of age. Dashed line marks 
the threshold of 100 Myr. Upper age limits are denoted by arrows whose hats correspond 
to the uncertainty in f^. When only an age range is known, no filled circle was plotted. 
Red symbols mark those debris systems whose existence was explicitely confirmed by an 
instrument independent of IRAS (Sect. 4.2). 
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Fig. 5. — Fractional luminosity of the infrared excess as a function of age. The shadowed 
bands mark the evolutionary models of Dominik & Decin (2003) (for detailed model param- 
eters see Sect. 4. 3). Upper age limits are denoted by arrows whose hats correspond to the 
uncertainty in f^. When only an age range is known, no filled circle was plotted. 
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Fig. 6. — Positions and proper motions of proposed new members of the Tucana-Horologium 
and GAYA2 associations (see Appendix B.) 



